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SUMMARY 

The paper presents the stability study of rotor/bearing systems. Even though it was 
limited to study of a fully lubricated bearing subject to oil whirl, and further 
limited to low eccentricity region for linearity and with only one type of lubricant, 
it can be seen that the perturbation methodology, together with the sorting of the 
impedance terms into direct and quadrature with respect to input force can be very 
useful to the general study of stability. Further, the concept of active feedback 
should assist to increase knowledge in rotor system stability. While there remains 
a large amount of study to be accomplished, perhaps some more tools are available 
to assist this field of analysis. 


INTRODUCTION 

The problem of rotor vibration in plain fluid bearing becomes more and more 
serious as operating speeds of turbomachinery continually increase. Two of the 
most significant sources of vibrations are fluid lubricated bearings and seals. 
Despite the high interest of field engineers and researchers which has followed 
the widespread application of oil /gas bearings and seals, the problems created by 
their hydrodynamic actions in turbomachines are still not solved completely 
[1-18]. 

This paper shows that the dynamic behavior of a rotor and its stability can be model- 
ed as an active feedback mechanism, which should greatly help in clarifying the 
stability problems of rotor system. It further shows testing and data techniques 
for bearings and seals, and the evaluation of the performance of a perturbation 
test bearing. 

The principal observed dynamic phenomena due to the lubricant action during shaft ro 
tation are known as oil whirl and oil whip. Oil whirl appears at rotational speeds 
lower than twice the first bending critical speed of the rotor (corresponding to its 
first natural frequency), as a forward circular motion of the journal with the 
frequency nearly equal to half the speed of rotation. This motion is related to a 
loss of stability of a pure rotational motion (about a central position of the 
journal in the bearing, or more generally - an equilibrium position, determined by 
the static load) and creation of a stable steady circular motion with an amplitude 
affected by nonlinear factors of the oil film. In the performance of real machines, 
oil whirl may cause serious machine damage since the level of the steady-state 
whirling amplitude may exceed admitted tolerances of vibration. 

Oil whip also has the character of forward motion, but unlike oil whirl, is locked 
to the frequency of the first self balance resonance (critical speed) of the rotor. 
The rotational speed must be of at least twice the resonance speed for oil whip 
instability to occur. Sometimes the rotor goes unstable in oil whirl, then converts 
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to oil whip when higher rotative speeds are reached. Oil whirl is not known to 
occur at higher speeds than oil whip. 

In earlier years it was believed that for all of the forward circular instabilities 
that the range of rotative speeds from 2 to 2.5 times the self-balance resonance 
speed was required. However, more recent works [1] show that such instabilities 
may occur at any rotative speed ratio with respect to balance resonances, depending 
on the type of instability. Some machines in the field with severe gas whip problems 
show that this instability cannot be cured with the machine rotational speed below 
the first self balance resonance. 

The literature on oil whirl/whip phenomena is quite rich [1-18]. The results con- 
cerning the loss of stability of the pure rotational rotor motion are usually pre- 
sented in the form of stability charts describing the behavior of the rotor in the 
vicinity of the steady equilibrium position for different variable parameters. The 
stability rules indicated by different authors rarely agree in detail. This is not 
surprising - rotor response is very sensitive to slight variations in parameters 
(e.g., a change of temperature by a few degrees, or minor geometric asymmetry in the 
bearing may stabilize the rotor), and it has not been recognized that the active 
feedback is involved. 

Classical analysis of the rotor/bearing systems usually starts from the very gen- 
eral Euler-Navier-Stokes equations. The solution of these equations, which may 
supply information concerning magnitude and functional relationship of hydrodynamic 
forces is unfortunately not easy to obtain. Analytical solutions require a lot of 
simplfying assumptions (such as idealized boundary conditions, neglecting thermal 
effects, ideal characteristics of lubricant, assumed pattern of the pressure distri- 
bution, neglecting curvature and capillarity effects of the lubricant, laminar flow, 
etc.). The resulting solution represents, then, only a rough approximation, as some 
of the neglected factors may happen to be dominant. 

Numerical solutions of the Euler-Navier-Stokes equations, requiring, just as much, 
a number of simplifying assumptions, involve a lot of time-consuming computations’ 

[6] and can solve only some particular problems. Progress in computer facilities 
and numerical techniques, however, improves the possibility of achieving solutions, 
which allow for rotor/bearing coupling of motion and further parametric analysis. 
Parallel to the progress in computer sciences, a remarkable progress in experimental 
techniques and measurement facilities has been observed. The experimental results 
bring better understanding of dominating phenomena occurring during rotor/bearing 
motion. This, in turn, will clarify eventually the simplifying assumptions intro- 
duced in theoretical analyses. 

In the present paper an experimental method of identifying the hydrodynamic forces 
in a plain cylindrical bearing during the oil whirl condition will be dicussed. 

The results are applicable to rotor/bearing or rotor/seal systems. The approach, 
based on analysis of forced response of the perturbed system, presents a continuation 
of previous research [7,8,9,12,13]. The mathematical model of the system considered 
is based on the linear theory of the rotor motion and fluid flow effects, as this 
series of tests have been deliberately run at low eccentricity numbers (oenerallv 
less than 0.5). 

The agreement of the experimental results and the predicted motion is very good. 

The method of investigation is shown to be very effective and sensitive to any 
parametric variation, so that it can be highly recommended in many applications. 
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TEST ROTOR SYSTEM AND THE PHYSICAL MODEL 
9 

The experimental rig to 
perform the perturbation test 
consists of two shafts with 
circular cross sections (Figs. 
1,2). The main shaft rotates 
at the constant speed lUp 
(being quite well below the 
first critical speed) in a 
rolling, pivoting bearing 
and in a plain cylindrical oil 
bearing with diameter D=2" 
axial length L=1.5" and 
radial clearance C=0.015". 

The second shaft rotating at 
perturbing speed uUp in a pi- 
voting bearing and a rolling 
antifriction bearing, is 
attached to the journal of 
the main shaft at the oil 
bearing. The trans- 
parent bearing, fed by blue 
dyed oil type ISO-VG-32 from 
an outboard tank with a 
constant pressure p=7100 
N/m^ permits direct obser- 
vation of the distribution 
of the lubricant within the 
bearing . 

Both shafts are driven by 
separate electric motors 
Figure 1. - Rotor test rig. 1 - Bearing, 2 - Main (q. 75 HP for the main rotor, 
shaft, 3 - Perturbing shaft, 4 - Perturbing q ^ perturbing 

speed controller, 5 - Additional springs, 6 - through a flexible rub- 

Main shaft tachometer, 7 - Thermometer, 8 - ber coupling. For compen- 

Proximitors. 9 - Perturbed Reservoir. sating weight of two shafts, 

a pulling vertical force is introduced in such a manner that during rotation of the 
main shaft, its response is a pure rotating motion in the center of the bearing. For 
this purpose, the main shaft is also straightened and balanced. At the same section 
of the shaft the additional symmetric springs can be introduced. The perturbing force 
is created by a controlled unbalance at the second shaft (calibrated unbalance 
masses attached to a single light disk mounted on the shaft). The unbalance inertia 
force is then proportional to the square of the perturbing speed utp . 

The motion of the shaft journal is observed by two displacement motion, noncontacting 
probes (vertically and horizontally located) plus a Keyphasor probe to provide ac- 
curate phase and speed signals. Additionally the constant rotational speed of the 
main shaft and slowly variable rotation speed of the perturbing shaft are recorded. 

A special device designed to control constant angular acceleration of the perturbing 
shaft was used. The Keyphasor, horizontal and vertical probe signals, filtered to 
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frequency, were continuously observed on the oscilloscope mon- 
itor (orbits and time variable displays) and stored in the computer. The computed 
then provides numerical results and graphical representation of the amplitude/phase/ 
frequency relationships, or their derivatives ("Direct" and "Quadrature" Impedances 
see expressions (9) and (10)). The physical model of the system is prLenS Ir ’ 


MATHEMATICAL MODEL OF THE SYSTEM 



Figure 2 . - Physical model of system. 


The mathematical model is 
based on the assumption of 
"small" deflections, allow- 
ing for linerization. The 
angular motions of the 
shafts are then approx- 
imated by their linear 
radial deflections. The 
hydrodynamic forces are* 
introduced in the oil 
bearing are introduced 
in their linear approxi- 


(x^’and vIrtTc^"(vrdffle"J "^' Vll"' of'f^'eedr(hoHzo^tal 

Lx; and vertical (y) deflections of the mam shaft, as seen by the probes) The 

thP^en supposed rigid. External damping is neglected. With th«e assumptions 
the equation of motion of the system is as followsf assumptions. 


where: 


Mz + (JiiUp + J2tUp)iz + Kz + F 
. 2 


^2 - 


= U e 
P 

U = (i|-Jl!2)J!. 
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z = X + iy, 
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i Pi , 
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r/£i£4 , 

1^4 
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(3) 

(y) deflections 
main (1) and 
, are cor- 



... „ va, lauic cue nunzoncai (x; and vertic 

of the journal; Ij, I I are moments of inertia of t 
perturbing (2) shafts (abidt thi x or y and z axis) £ . . . . 
responding lengths (Fig. 2); u)„ is main shaft rotational spe 
shaft rotation speed, K - stiffness coefficient of the addit;^^a?spr1n«-'m''''' 

fo^cr r 1 m dhbaiance, P - vertical constant polling"’ 

force, mi. ">2 - masses of the main and perturbing rotors, g- gravity acceleration 
F represents the hydrodynamic force in the bearing: gravity acceleration. 

F = (m r n.ji)z r (d^-dji)i + 

where and are bearing coefficients. 

lake^lIto'T^s^derlMor Mrodynamic forces is 

SOLUTION OF THE MATHEMATICAL MODEL 

We simplify the mathematical model (l)+(4) neqlectino ovrosrnmV to.^mc ri -i -n^ 
and the coupling "fluid Inertia" tern, (1^4) ® avro.copic terms (J2-J3-O) 
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If the condition of stability is satisfied, 

where f, f^ are natural frequencies of the free vibration (see Appendix for 
detail s)i then the dominating solution of (1) will be the forced response: 


1 (U) t+y) 

z = A e P + B , 

(6) 

where amplitude A and phase angle y are 


2 

A = U cOp /(D^ + 

(7) 

Y = arctan (- Q/D) , 

(8) 

with 


Direct Impedance: D = -(M + m) uj ^ +K 

(9) 

and 

Quadrature Impedance: Q= d^uip - k^Wj^. 

(10) 



Note that the phase angle y 
corresponds to the attitude 
angle of the bearing. 

The static displacement in the 
solution (6) is 


figure 3 . - 
responses 


Phase frequency 
of system. 


and amplitude/frequency 


B = 




( 11 ) 


— ? 

(K+k^l'^+k'^ 




It will serve for evaluation of 
the stiffness coefficients by a 
preload test with no perturbation 
(u) =0). During the perturbation 
te?t (variable lo ) the pre- 
load was eliminated (Pi=0). 

TEST RESULTS. IDENTIFICATION 
OF THE IMPEDANCES AS FUNCTIONS 
OF ROTATIONAL SPEED. STABILITY 
CHECK 

The experimental tests produced a 
set of results for the response 
amplitude A and phase angle y 
versus (slowly varying) per- 
turbation speed. The para- 
meters of the set were rota- 
tional speed and temperature. 

The tests were performed for 
three cases: with additional 


springs (System "S"), with additional springs and mass (System "SM") and without 
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them (original system "0"). Fig. 3 presents an example of the test results for 




REVERSE 


ORIGINAL SYSTEM (0) 

WITH ADDITIONAL 

SPINGS (SI 

WITH SPRINGS 
& MASS ISM) 




PERTURBATION SPEED co. 
(rad's 100] 


^<5_H(0) 3 

I 

rH{Sl ; 

I 


FORWARD 


H HORIZONTAL 
V - VERTICAL 


H(SM)- 


2jo>« 3 

PERTURBATION SPEED u;„ 
[rad/s 1001 


rf = 64 4c poise (T - 67 '-’Fl 
COr - 200 rad/s 


Figure 4. - Impedances versus perturbation speed 
for three cases . 


a rotational speed u)„ =200 rad/s. 

For a perturbation speed slightly 
less than half of the rotational 
speed, a well pronounced resonance 
was observed. This pattern repeated 
for all values of the rotational 
speed. In the resonance vicinity, 
the phase angle sharply changes 
values between zero and -90°. 

This behavior corresponds well to 
the predicted model of the sys- 
tem (see expressions (7), (8)). 

The system without additional 
springs and mass shows the highest 
resonance amplitude. An addition 
of springs lowers the amplitude 
considerably. 

For the reverse perturbation speed 
resonance has not been observed 
although the system attempts 
to whirl forward. The ampli- 
tudes of vibration slowly in- 
crease with negative frequency, 
the phase angle slightly changes 
at the frequency corresponding 
to minus half u>n. Major differ- 
ences in horizontal and vertical 


tnree cases. phase angles are due to the 90° 

fSequeL^Tegion. nonsymmetricity in the system, especi^ly in^'tL'"-!/^^^"' 

IHd Q^adr'iture'lmped‘anL‘sr''‘ Direct 


D = U oj COSY /A , 
P 


Q =-U (j siny /A 
P 


generated relationship between Impedances and perturbation 

D?nlteS iS “■ ceerespond to the results 

J ?• ^ ^ ^ symmetry of the system for the forward 

perturbation speed was maintained - the amplitudes and phases showed to be the same 
for vertical and horizontal responses, the orbits were circular. At the reverse 
perturbation speed, we observed differences in horizontal and vertical Impedances. 

Figure 5 presents a set of results of the response of the original system "0" for 
variable rotational speed. All other parameters of the system were keot constant 
during this test. The Direct Impedances have a roughly parabolic shape while the 

Quadrature Impedances are nearly straight lines. uoi.c snape, wni le the 
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The data taken from Fig. 5 served to 
obtain the cross relation for the 
Direct Impedance versus rotational 
speed with the perturbation speed 
as a parameter (Fig. 6). An analy- 
tical approximation to the Direct 
Impedance as straight line function 
of the rotational speed has been 
adapted: 

D — C jUJp- C 2 • 

The next step of identification was 
an evaluation of the coefficients C, 
and Cp as functions of the perturbation 
speed w . For this purpose, the val- 
ues of Ehe Direct Impedance at zero 
perturbation speed and the slopes 
of the straight lines present in Fig. 

6 were crossplotted as functions of 
the perturbation speed. 

For both of these curves the parabolic 
approximation has been adapted [15]: 

Cl = - 1.8-10~^ u) 2 - 2.5 u) -200 [kg/s] 

Co = + 0.25 u) 2-950 u) ~ 9 . 5 - 10 '^ [N/m] 
P P 


Figure 5. - Impedances versus perturbation 
speed for variable rotational speed. 


D/0.77 = -(1.8-10'^ u)d + 0.25)w 2 - y, 

K p p 


Finally the Direct Impedance is describ- 
ed by for the following analytical 
form: 


(2.5wp-950)-200 + 9.5104, (13) 


where 0.77 was not yet taken into account the dimensional factor. Comparing this 
expression with the expression (9) derived from the mathematical model of this 
system and introducing M=0.4 kg, K=0, the identification of the parameters of 
the rotor can be completed: 


m = m(u)j^) = (1.410 - 0.21) [kg] 

dt= d^(u)j^) =-(1.9 - 730) [kg/s] (14) 

k = k (u)„) = (7.3-10^ - 150 u).) [N/m]. 
r r K K 


The Quadrature Impedances (Fig. 5) were identified as the straight lines; varying 
with rotational speed. 


where 


Q = (970-0. 69wp)(u)p-Au.^), 

“ p1q=o 


(15) 

(16) 
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is the ratio of the pertur- 
bation speed for which the 
Quadrature Impedance has a 
zero value to the rotational 
speed u)„. The parameter \ 
appeared to be nearly constant 
for all performed tests: 

A = 0.48. 

Its value agrees with the oil 
whirl ratio observed in prac- 
tice. The value Aujp corre- 
sponds to the average oil 
velocity. 

The comparison between the 
expressions (10) and (15) 

R r >4 4 . • .i leads to the identification 

Figure 6. - Direct impedance versus rotational speed of the system parameters: 
and variable perturbation speed. 

d^ = d^(io^) = (970-0.69 u)j^) [kg/s] 



and 


*^t “ ~ ~ ***R^ [kg/s] = Ad^(u)p) 

"®t ~ has been already assumed. 


(17) 


As we notice, all parameters of the system are functions of the rotational speed uj 
It was not possible to obtain test results for very low values of the perturbation^ 
speed. For this reason a "displacement" test was performed. The main shaft rotating 
at constant speed was perturbed by a constant pulling force P. The probes measured 
vertical and horizontal displacements. The force-response characteristics showed 
to be linear in relatively wide range of displacements, however, their slopes varied 
considerably with the rotational speed. The formula (11) served for the identi- 
fication the parameters k^ and k^ (K=0 in this test*): 

k = 4000 N/m (for 0 < u)„ ^ 150 rad/s) ri8l 

k[ =( 490-0. 27uj^) [kg/s]. ^ 

The evaluation of the coefficient k^ was not very precise. However, its relatively 


The calculation of the coefficients k and k. was performed indirectly. First 
the slopes of the straight lines correspondingHo the pulling force P, versus 
horizontal (x) and vertical (y) displacements were evaluated for different 
values of the rotational speed. 

The slopes are equal to "directional" stiffnesses K =P^/x, K =P,/y, corresponding 
to X and y. Taking (6) and (11) into account the stiffness ^^oerf icients 

k = K / (1+(K /K )^), k^(o„ = k K /K 

r y yx tR ry'x 

were then calculated for each value of the rotational speed. 
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small value was noticed especially for higher lUp . The agreement with the previously 
evaluated coefficient (expression (17)) is very good. 

Having evaluated the parameters of the system the stability check has been 
possible. 

Introducing the expressions (14) and (17) into the stability criterion (5) and 
rearranging terms, we obtain the cubic polynomical inequality: 

1.29 + 3030 . 6.17-10^u)^-2.29 10® < 0 (19) 

The polynominal (19) has one positive root lUp = 366 rad/s. The stability 

During the experiment it 
was observed that for rot- 
ational speeds higher than 350 
rad/s the rotor became un- 
stable, which is consistent 
with the above analysis. 

We notice that the value 
of amplitude peaks is close- 
ly related to the difference 
between perturbation speeds at 
which Direct Impedance and 
Quadrature Impedance equal 
to zero (Figs. 3 and 4). 

This difference is referred to 
as the stability safety margin 
A (see Appendix). In Fig. 7 
ROTATIONAL SPEED wnirad/si the values of perturbation 

speed for zero impedances 

Figure 7. - Perturbation speed at zero impedances and been plotted. The mar- 

stability margins. gin of stability decreases 

with the rotational speed. For the considered systems ("0","S", and SM ) the 
zeros of the Quadrature Impedance remain the same. The additional springs and 
mass have influence on the Direct Impedance only. 

At the reverse perturbation speed, the Quadrature Impedance is negative for 
all speeds (compare (10) for u) negative, see also Figs. 4 and 5). It does 
not cross zero; the system is very stable. 

Figure 8 presents the impedances for the squeeze film test («>n=0). We notice the 
symmetric parabolic shape of the Direct Impedance and antisymmetric (straight line) 
for the Quadrature Impedance. (Slight displacement of the lines from zero at 
u) =0 is due to very low input force level). As expected, tangential stiffness, 
bBing related to the rotational speed, equals zero when ujj^=0. 

In the squeeze film test the Quadrature Impedance contain then only the radial 
damping term d (see (10)), which appears to be positive. The condition of 
stability (5) ’fs satisfied for all rotational speeds. The displacement of the 
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-horizontal 


- 1 000 

inch 

-2 


n ^ 64.4c poise (T ^ 67 °F) 
= 0 (Squeeze Film) 


■ ^ Direct Impedance parabole from 

the symmetric position for iUp=0 

- — I I I ^ o ] 2 3 into the region of positive ^ 

pJoRBAxIoNspLci: perturbation speeds for u, >0 

vert' I indicates a relationship of the 

N. tangential damping d with rota- 

// S 5--2 within the 

' 1 investigated range of speed. 

a (See (14) for comparison). 

{C --3 horizontal^ \ ' 

Q 

n 64 4c poise IT = 67 °F) The Fig. 9 presents the Direct 

I ^-^=j^oee.eF,M ^ gpd Quadrature Impedances for 

variable oil temperature. The 

variations of the Quadrature 

Impedance are quite significant 
O ^ horizontal-^^^ while the Direct Impedance 

PP.,PP I “- 2 5 remains nearly invariable in 

tKbt - 1000 -^ vertical temperature range, espe- 

3 -2 1 o 0 cially for lower values of the 

' f ' + ^ ^ \ 1 1 , [ ^1 perturbation speed. Let us 

perturbation SPEED j evaluate the Quadrature Impe- 

I, ad/s 1001 ’ dance as a function of tem- 

--25 perature. From the graphs Fio. 9 

FORWARD the relation between the slope 

of the Quadrature Impedance 

' iso J versus temperature and the 

FiQure 8 - Tmnactan/.... speed ratio \ can be obtained, 

squeeze* film ^ versus rotational speed for Finally the Quadrature Impedance 

n/n - 7 ^ evaluated as follows [15]: 

Q/0.77 = (1700 - lOT) [u)p-A(T)u,^J [N/m], 

A(T) = 0.48 (1.1 - 1.8 lO'V), 

where T is temperature in Farenheit degrees. Further we derive the relationshin 

s^that th^Srl? /■'scosity (measured by Cannon-Fenske Viscosimeter^OO) 

?? Quadrature Impedance can be expressed in terms of the oil viscosity 

he Direct Impedance is known for any temperature, then a correcting factor 

temperaJuJe°U5L Impedance corresponding to any different 

Summarizing the test results, we see that the predicted model well reflects the 
features of the system. The Direct and Quadrature Impedances can be 

ACTIVE SERVOMECHANISM CONCEPT 

a lubricated bearing to support both steady-state and dynamic loads 

well known that'll consisting of a fluid wedge is well known. It is also 
ell known that a 360 lubricated bearing has a propensity for instability known as 

no: "'ll- ”“‘=‘="‘''''9 work of Newkirk and Taylor 
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Figure 9. - Impedances versus perturbation 
speed with variable viscosity (temperature) 
of the oil. 


The Study of rotor instability has 
generated a vocabulary including 
such expressions as cross springs, 
cross damper, (and even cross masses), 
negative damping, etc. Such cross 
actions exist as effects, unlike gyro- 
scopic action which is a genuine cross 
coupling physical law. The reason 
why these cross coupling bearing and 
seal effects occur is that such ele- 
ments clearly have a behavior much 
more complex than any possible passive 
systems. These machine elements 
contain active feedback, in which 
exists a mechanism transferring rota- 
tional energy into lateral vibration 
energy. This is true for every identi- 
fiable rotor instability of the forward 
circular whirl and whip category, 
including shaft internal friction, 
aerodynamic cross coupling or steam 
whip. In particular, the behavior 
of a 360° lubricated bearing is 
analyzed herein. The situation 
in any gas or liquid seal is also 
generally in the 360 degree category 
so that the results can then be ex- 
tended to any seal or bearing dynamic 
behavior. To show in simple descrip- 
tion that such systems are actually 
closed loop servomechanisms, observe 
the performance of the bearing in 
response to a steady- state load. 


Apply a downward pressure to a shaft rotating anticlockwise with no initial 
dynamic or steady-state preload. The obvious and well known result is that the 
shaft moves downward very slightly, and moves strongly to the right. This 
motion, of course, provides the converging lubricant wedge to support the 
applied downward load. Since lubricant has virtually no spring effect of its 
own, the action generates the spring cushion synthetically. Applying a down- 
ward steady-state load as a step function, the first effect is for the shaft 
to move downward in response to passive direct (radial) damping (d ) • The 
result of this downward action is to create a restriction of oil flow at the 
bottom. This in turn creates a converging wedge on the left, which forces the 
shaft to the right. This process continues (usually with some damped whirling), 
until the shaft reaches an equilibrium position in which forces are balanced. 
This position to the right is sufficient to create the constriction which 
provides the wedge to support the shaft against the input force because the 
input to the active element is motion and the output is force, as previously 
described. 


The dimensions of the active elements are that of impedance. It then seems 
reasonable to present the rotor/bearing system as a servomechanism with a feedback 
loop containing active elements of the system (Fig. 10). The general transfer 
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function of the system (described by Eqs. (1) and (4)) is given in the Appendix. As 
it has been shown in the previous section, the parameters related to the bearing 
depend directly on rotational speed - they represent active elements. The term 
"active" refers to the mechanism of transformation of rotational energy into vibra- 
tion energy, which provides the major reaction to input forces. 

The input to the system is represented by the exciting periodic force as well as 
static load ((|>(iUpi) - see Appendix). 
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Figure 10 


The output is the motion - 
vertical and horizontal dis- 
placements of the main shaft 
(z(iu)p)). 

The output motion with respect 
to the input force shows all of 
the active and passive terms in 
simple series (see (A8)). As 
such, it could be optional 
which ones are put in the 
active category, with the 
remaining ones left in the 
passive category. However, 
the static loading test pro- 
vides a perturbation speed of 
ui_=0 so that all terms con- 


- Rotor- bear i ng active feedback diagram. 

taining to as a multiplier are deleted. These terms deleted entirely are the 
passive terms. Further, the squeeze film tests, where rotative speed tUp=0 deletes 
all terms containing lUp as a multiplier. The terms which are deleted entirely can 
considered as the active terms. These observations indicate the proper location of 
the terms. 



Figure 11. - Force diagram for 

rotational speed below resonance. 

one moving part, as well as extremely 
performs several functions. 


If a further check is desired, since it is 
fundamental that springs with negative 
stiffness cannot exist as real passive 
impedances, any negative stiffness term 
component necessarily represents the active 
portion of the feedback loop. Since spring 
force is motion times spring coefficient, 
in order to supply the restoring force 
against input force the elastic force must 
be at nearly 90 degrees with respect to 
motion, because the motion itself is nearly 
90 degrees with respect to input force 
(Fig. 11). 

As a result, the so called "cross spring" 
is actually the very real (but syn- 
thetically manufactured by feedback) 
restoring spring force. In review of the 
active behavior of such a bearing it is at 
the same time extremely simple with only 
complex, because this one working part 
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(20) 


Summarizing these considerations we describe the transfer function of the active 
feedback loop in the following form (Fig. 10): 

^ — j OC 

k^-iu)pk^ sirvx =k^iUp(cos « -i sin « )=k^u)j^e , 

where k^=k^cos«, k^ =k^ + k^ , sin « =A.d^, 

The active loop then contains a fraction of the u)„ dependent coefficient k^ 

(radial stiffness) and full tangential stiffness coefficient. The last one is 
entirely the result of oil wedge formation. As it was identified in the previous 
section, the radial stiffness contains a constant term and a term varying with the 
rotational speed. However, a precise identification was difficult in view of its 
relatively small value. 


It is supposed that the terms located in the feedback loop are directly responsible 
for the transfer of rotational energy into vibrations. The angle « represents the 
angle between the observed motion and the active force (Fig. 11). This angle is 
usually a little less than 90 degrees for 360® fluid bearing. For other types of 
bearings it is nearer 0 degrees (e.g. for sleeve bearings with stable void islands in 
the thick film region, generally referenced to as "180 Degree or Half Sommerfeld 
bearings). 


In the present study it has been found that the radial spring coefficient has a low 
value, for « =90° cos « is very small - this is reflected in the expression (20). 
However, at this time we have not found a means to measure directly the angle «. 


All other terms of the rotor system are located in the "passive" part of the 
servomechanism model (Fig. 10). Being functions of the rotational speed, they may 
also modify the rotor response in some "active" way. however to a lower degree than 
the ones related to stiffnesses. 


A rotating shaft in a bearing does (1) have passive impedance, (2) transfers energy 
from input torque by dragging lubricant around. (3) moves sidewise of the input 
force (as best it can) to provide a restriction of flow to the lubricant to build 
a converging wedge for support, and finally (4) synthetically builds a force to 
counteract the input force. 

Rotor systems, as currently configured, are rarely well stabilized because (1) their 
basic algorithims are complex, nonlinear and nonsymmetric, and (2) there has been 
lack of recognition and knowledge of this feedback nature of rotor behavior. The 
present paper throws some light on these problems, however much more research is 
needed to provide answers to all questions related to the rotor/bearing/seal 
stability. 
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APPENDIX 

AI. EIGENVALUE PROBLEM FOR THE EQUATIONS (1) AND (4) 

The homogeneous equation (1) with hydrodynamic forces (4) has the following 
eigenvalues: ^ 

^I’2’3’4 ^ ^r* i (d^-JjWp-J2Wp) )/(2M+2m) + (aj^+ ia 2 )j 

(Al) 

where 

a- = [(-1)’^ ((d^-JjOOp-J2(iJp)^ ■ ^r 4( (M+m)k^+m^k^a)p))+ 

/ 2 2 7 

+ C( “ dy^ + 4( (M+m) ky'+m-j-k-f^cop) ) + 

' / A O \ 


Oil 

.4 (d^(d^-J^a)p-J2Wp) + 2(m^k^-(M+m)k^a)^))7"| "V(v^'(M+m) ) , j=l,2. 


The condition of stability (negative real parts of the equinvalues (Al) i 
follows: 


s as 
2 


d^{(M+m)k7m^k^Wp)-d^{d^-J^c0p-J2t0p)(m^k^-(M+m)k^Wp)-(m^k^-(.M+m)k^a)p) > 0. (^3) 


For Jj - J 2 = m^ = 0 it reduces to 
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(A4) 


-(M+m) (k 


t'^R 



or after transformation to 
-Ifal < < Ifil , 

where 




d^/(2M+2m) +(-l)'^((d^/(2M+2m)) +k^/(M+m))] 


(A5) 


, j=l,2 (A6) 


are natural frequencies of the system at the limit of stability, i.e. , when 
instead of inequality (A4) we have the equation. As we can easily check fi and f 2 
are also the roots of the Direct Impedance (9). 


For k^ > 0 the quantity 

A = fi - u)pk^/d^ > 0 

gives a margin of safety for the stability. As the value Wpk^/d^ is the root 

of the Quadrature Impedance (10), the stability of the perturbed system is 
assured if the zero of the Quadrature Impedance appears at a lower perturbation 
speed ujp than the zero of the Direct Impedance. 

All. FOURIER TRANSFORMATION OF THE EQUATIONS (1) and (4). 


The mathematical model of the system can be presented applying the Fourier 
transformation and the control theory formalism. From (1) and (4) we get 
then 

z (iu)p) [D (iujp) + iQ (i^p)] = L [ 0 (iWp)] 

where z (i“Jp) is Fourier transformation of the variable z(t) (2), iiUp is 
the Fourier operator, D(iujp), Q(iujp) are given by (9), (10), L [ 0 (iWp)] is 
the Fourier transformation of the input forces 0 = Uio^ e - iPi- [D(iu)p) + 

+iQ (iw)p)] is the transfer function of the system. 

NOTATION 

A Amplitude of the response _4 

C Bearing radial clearance (0.015"=3.8 10 m) 

d ,d Radial and tangential bearing damping coefficients 

^ Direct Impedance 

g Acceleration of gravity 

I ,1- Moments of inertia of the main (1) and perturbing (2) shafts 

^ about X or y axis (24.8 gm^ , B.Tgm^; 1^=50 gm^ for the case of 

additional mass) 
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Moments of inertia of the main and perturbing (2) shafts about 
the axis of rotation (O.OSSOgm^.O.OSgm^; I, =0.0652gm2 for 
case of additional mass) 

Radial and tangential bearing stiffness coefficients 
Length of the perturbing shaft (9.5"=0.241m) 

Distance from the center of bearing to the probe (0. 625"=0. 0159m) 

?Innfh‘'«/rr • to the pulling force section (2"=0 0508m) 

Length of the main shaft (11. 75"=0. 298m) u.uouom; 

Distance to the position of the perturbing force (7. 7"=0. 195m) 

asses 0 the main and perturbing shafts with accessories (280a 1500 ") 
Perturbing unbalance mass (4g) les uoug, ibOg) 

Radial and tangential bearing "fluid inertia" coefficients 
Original system 
Pulling force 
Quadrature Impedance 

Radius of the controlled unbalance (33/32"=0. 0262m) 

Radius of the shaft journal (1. 00"=0. 0254m) 

System with additional springs 
System with additional springs and mass 
Phase of the response - Attitude angle 
Stability margin 


Dynamic viscosity of oil (poise) 

Oil rotation speed ratio 
Perturbation speed (perturbing shaft) 
Rotation speed (main shaft) 


322 



